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We examine the problem of constructing spectral representations for two point correlation func- 
tions, needed to write down the QCD sum rules in the medium. We suggest constructing them from 
the Feynman diagrams for the correlation functions. As an example we use this procedure to write 
the QCD sum rules for the nucleon current at finite temperature. 
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Q ' I. INTRODUCTION 

<n : 

(3JQ| It was realized quite early in their formulation that the QCD sum rules in medium have several new aspects 
not shared by the corresponding vacuum sum rules H . On the operator side, there appears in general new operators 
' in the short distance expansion of the correlation functions, due to the presence of the four-velocity vector of the 
medium ^,|^]. Even for the old operators, the evaluation of their matrix elements in medium may need further 
untested assumptions. As an example, the four-quark matrix element in nuclear medium is evaluated by the ground 
state saturation, different from the vacuum saturation used for the vacuum matrix element Again on the spectral 
H ' side of the sum rules, unlike the vacuum case, a communicating single particle (with medium dependent mass and 
coupling) alone does not saturate the spectral function in the low energy region; in the medium relevant two-particle 
04 i states are equally important. [Throughout this work we shall be interested only in terms linear in the distribution 
CO ■ function. For higher order terms, contributions from multi-particle states need be included.] 

The present work concerns the spectral side of the QCD sum rules. We take a closer look at the contributions of 
the two-particle intermediate state. If we draw all the 1-loop Feynman diagrams for the correlation function, we see 
that in addition to the diagrams with the (single particle) pole alone and the (two-particle) branch points alone, there 
are other (1-particle reducible) diagrams, which appear as the product of the pole and the branch points. They have 
remainders after extracting contributions which modify the pole parameters. It is these remainders which are not 
always included in the saturation scheme, even though they are generally of the same order as those it retains. 
We illustrate our point by discussing the nucleon QCD sum rules at finite temperature T. This choice is dictated 
Oh, by the simplicity of the nucleon current correlation function at finite T and the availability of results for comparison. 
For T < fx, the pion mass and zero nucleon chemical potential, pions dominate the heat bath. The leading term in the 
effective chiral Lagrangian involving nucleon and pion fields describes the interaction well in this case, allowing the 
calculation of shifts in the pole position and the residue [Q . The operator side of the sum rules is also simpler at finite 
T than say, in the nuclear medium: To 0{T 2 ) only the Lorentz scalar operators entering the vacuum sum rules can 
contribute, the contributions of the new, non-scalar operators being at least of 0(T 3 ). Further, the thermal average 
of the operators can all be reduced to their vacuum expectation values, so that no further assumptions are required to 
evaluate them, other than those needed in the vacuum case. What is more, these sum rules have already been worked 
out H and it would be interesting to see if the above-mentioned remainder terms bring in any new contribution. 

For increased sensitivity, we work with the so-called subtracted sum rules, i.e. the finite T sum rules from which 
the corresponding vacuum sum rules have been subtracted out, displaying only the contributions of 0(T 2 ). The 
spectral representations for the sum rules are worked out for physical masses and couplings. Only in the evaluation 
of the resulting integrals, do we take the chiral limit. For completeness we have derived in an Appendix the Wilson 
coefficients for the product of nucleon currents using a simple coordinate space method. 

In Sec. II we describe the construction of the spectral representation for the correlation function. In Sec. Ill we 
evaluate these contributions in the chiral limit. Sec IV displays the sum rules and their evaluation. Sec.V contains 
our concluding remarks. Appendices A and B respectively give the absorptive part of 1-loop diagrams at finite T and 
the Wilson coefficients in the short distance expansion of two nucleon currents, along with the thermal averages of 
the relevant operators. 



X 

5-H 



*Email address: mallik@tnp.saha.ernet.in 
^Email address: sourav@veccal.ernet.in 



1 



II. CONSTRUCTION OF SPECTRAL REPRESENTATION 



Consider the two point correlation function 

H(p) = i / d 4 xe lpx Tr{pT V (x)f](0)) (2.1) 



of the nucleon current rj{x)D,i, composed of three quark fields, having the quantum numbers of the nucleon, the indices 
D, i referring to spin and isospin [p|flO|]. Here p = e~@ H /Tre~@ H is the thermal density matrix with the Hamiltonian 
H of QCD. For the spectral representation we do not need to spell out the quark structure of r](x). 

A formal spectral representation for the two point function in E = po at fixed p may be obtained immediately 
|TT| , fl2"f . Evaluating the trace over a complete set of eigenstates of four- momentum and then again inserting the same 
states between the two currents to extract the x-dependence, one arrives at 

U(E, p) = iImIL n {E, P1 + — J E' - E ' ^ ' 

where the resulting double sum over states within the integral may be converted back to the Tr(ace) to get, 

M u (£,p) = \ J cPxe^Tripfrix), fj(0)]) (2.3) 

Here the index 11 in I mil reminds us that we are working in the real time formulation of the finite temperature field 
theory, where all two point functions assume the form of a 2 ® 2 matrix. Thus the expression (2.1) is actually the 
11-component of the corresponding matrix. However, no information is lost by considering the 11-component only, as 
for p 2 (= E 2 — p 2 ) space-like, where one writes the QCD sum rules, there are no other independent components. We 
also note that the imaginary part of the 11- component may be written as tanh(j3E /2) times a function, the latter 
coinciding with the imaginary part calculated in the imaginary time formulation 113] . It is thus convenient to define, 

ImU(E,p) = n- 1 coth{PE/2)ImU 11 (E,p). (2.4) 

We shall write the sum rules for p = and in the rest frame of the heat bath. Then we can decompose the two 
point function as 

U(E,p = 0) =Tl 1 (E 2 )+ lQ EIl 2 (E 2 ), (2.5) 
where the scalar functions have, in the notation of Eq.(2.4), the spectral representations, 

n,<^) = /!~^>, i = 1 ,2. (2.6) 

the integral running over two branch cuts, the short cut, < E 2 < (m — p) 2 and the unitarity cut, E 2 > (m + p) 2 , 
where m and p are the masses of nucleon and pion. 

Eq.(2.3) does not prove convenient to calculate the spectral function. Instead, we wish to identify the set of 
relevant Feynman diagrams and calculate the spectral function directly from them. In their work on the propagation 
of a nucleon through a heat bath, Leutwyler and Smilga [(?J considered the set of 1-loop Feynman diagrams for the 
correlation function, correcting the vacuum amplitude, shown in Figs. 1-4. From these diagrams they calculated the 
amplitude in the vicinity of the nucleon pole to obtain the shifts in the nucleon mass and the residue to 0(T 2 ). The 
same set of diagrams also suffice to saturate the spectral functions to the same order in T in the low energy region. 
Once the spectral representations of these diagrams are constructed, one can, of course, immediately read off the 
position and residue of the nucleon pole. 

It is simple to find the spectral representations given by these diagrams, except for a complication with diagrams 
of Fig. 1(b) and Fig. 2(a) and (b). The latter diagrams have the structure of the product of a (simple or double) pole 
times branch cuts, due to N and ttN intermediate states. We have to write these contributions as a sum of the pole 
and the branch cuts. In the rest of this section we perform this separation. 

Let us begin with the diagrams of Fig. 1 . Omitting the r)N and fj N coupling constants (to be reinstated in the next 
section), their sum is 

' ' -AE)-^—^---, (2.7) 



jqE — m jqE — m jqE — m 
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where the dots indicate the inclusion of the series of 1-particle reducible self-energy diagrams. The self energy matrix 
itself may be decomposed as 



£(£) = A 1 (E 2 ) + l0 EA 2 (E 2 ). 



(2i 



Given the absorptive parts, to which we turn in the next section, the scalar functions Ai(E 2 ) have spectral represen- 
tations like those for IT in Eq.(2.6). 
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FIG. 1. Nucleon pole and self-energy insertion. 
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FIG. 2. Vertex correction from ttN intermediate state. 
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FIG. 3. nN intermediate state. 



FIG. 4. Constant vertex correction. 



The second term in expression (2.7) is a product of a double pole and cuts contained in S. To write it as a sum of 
the pole and the cut contributions, we write T,(E) as 

£(£) = A 3 (E 2 ) + ( l0 E-m)A 2 (E 2 ), A 3 (E 2 ) = A^E 2 ) + mA 2 (E 2 ), (2.9) 

and expand the scalar functions around E 2 = m 2 to the required order to get 

E(E) = a+(E 2 - m 2 )b + ( 7o £ - m)c + R{E), (2.10) 

where, for short, 

a = A 3 (m 2 ), b = A' 3 (m 2 ), c = A 2 (m 2 ), (2.11) 

the prime denoting differentiation with respect to E 2 . The spectral representation for the remainder R may be 
obtained from those of Ai(E 2 ) by equating the two expressions (2.9) and (2.10) for S. We get 
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R(E) = (E 2 - m 2 ) 2 A 3 {E 2 ) + ( l0 E - m)(E 2 - m 2 )A 2 (E 2 ), (2.12) 

where 

x (F 2, f dE' 2 ImA 3 (E' 2 ) 

ME ) - J {E , 2 _ m2)2[E , 2 _ E2) (2-13) 
r , 9x f dE' 2 ImA 2 (E' 2 ) 

Inserting Eq. (2.10) in Eq.(2.7) and summing the series after isolating the contribution from R, we get the desired 
spectral representation for the diagrams of Fig.l as 

" 7 oJ-7m+ C a ) ~ m2)ME2) (7oi? + m )i 2m ME 2 ) + A 2 (E 2 )} (2.15) 

Next consider the diagrams of Fig. 2, each of which is a product of a simple pole and cut contributions. They are 
of the form 

- HE), (2.16) 

■foE — m 

where A(E) is the vertex with irN intermediate state. Again we decompose A(E) as 

A(E) = B^E 2 ) + lQ E B 2 (E 2 ) 

= B 3 (E 2 ) + ( l0 E-m)B 2 (E 2 ), B 3 (E 2 ) = B^E 2 ) + mB 2 (E 2 ), (2.17) 

and expand the scalar function B 3 (E 2 ) around E 2 = m 2 to get 

A(E) =d + S(E), (2.18) 

where d = B 3 (m 2 ) and the remainder S is 

S(E) = ( l0 E - m){B 2 (E 2 ) + ( lo E + m)B 3 (E 2 )} (2.19) 

where 

T'2 7™ D / I?'2 



B 3 (E 2 ) 



f dE' 2 ImB 3 (E' 2 ) 
J (E' 2 - m 2 )(E' 2 - E 2 ) [ ' 



The expression (2.16) now separates into the pole and the cut contributions 

d 



7o-E — m 



B 2 (E 2 ) - (-f E + m)B 3 (E z ) (2.21) 



Clearly to find the nucleon pole position and the residue, the terms with Ai(E 2 ) and Bi(E 2 ) in Eqs.(2.15, 2.21) 
are not necessary. But to write the QCD sum rules, which require the correlation function for all (large, space-like) 
values of E 2 , the entire spectral representation is needed. 



III. EVALUATION OF SPECTRAL REPRESENTATION 



Having separated explicitly the contributions of the TV-pole and the ttN cut arising from diagrams of Fig. 1-2, we 
proceed to calculate the spectral functions of all the diagrams. We need to know the interaction vertices present in 
the diagrams. The irN interaction Lagrangian is given, to leading order in chiral perturbation theory, by the familiar 
term, 

where qa = 1-26, the axial vector coupling constant of the nucleon and = 93 MeV, the pion decay constant. 
The coupling of r\ (x) with nucleon is defined as 
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(0\ Vl (x)\N 3 (p)) = 8 l0 \u{p)e-^ x (3.1) 

Its couplings with nucleoli and any number of pions introduce no new constants. They may all be reduced to the 
coupling with nucleon alone by using PCAC and current algebra 0). We conveniently represent these couplings by 
writing rj(x) in terms of the relevant physical fields, 

T)(x) = Xip(x) + ^-r • $~f 5 ip(x) - -^4>-$4>(x). (3.2) 

We now begin with the evaluation of the self-energy integral £(p). It is of the form of Eq. (A.l) of Appendix A 
with f(p, k) given by f(p, k) = ci^^(f> — ^ + m)^7s, c\ = — (3/ A)(gA/ F^) 2 . Note that any term ~ 7 • k in / will 
be zero after integration over the angle of k. Then for p — and in the chiral limit, it simplifies to f — 2c\k 2 .E^Q. 
Referring to Eqs. (2.8, A. 9) we get 

ImAi = 0, ImA 2 = — -y=7ii(|wi|)|o;i| , u x = — — — 

on both the branch cuts, < E 2 < m 2 and E 2 > m 2 in the chiral limit. The constant 2mb is given by 

f dE 2 ImA 2 (E 2 ) g 2 A T 2 . . 

The other two constants a and c are each of 0(T 3 ). 

The integrals in the remainder R may be evaluated in the same way. Consider As{E 2 ) given by Eq.(2.13). Antic- 
ipating the QCD sum rules, we actually need the Borel transform of the spectral representation at large space-like 
values of E 2 . It is given by 

„ T /w2 „ 2m 2 fdE 2 e- E2 / M2 ImA 2 (E 2 ) 

2m A 3 (M 2 )\ Borei - 1 



AP J {E 2 -m 2 ) 

= 9A_ e -m 2 /M 2 ll (3 4) 

16M 2 F*' 1 ' 

where M is the Borel mass. The Borel transformed integrals for (E 2 — m 2 )A^(E 2 ) and A2(E 2 ) are each of 0(T 3 ). 
The vertex A(E) is again given by Eq.(A.l) with 

f a (p, k) = c 2 ^7s(^ - ft + m)75, fb(p, k) = c 2 7 5 (^ - ft + m)^ 5 , c 2 = -(3/4)g j4 (A/i ;i 7r ) 2 , 

for diagrams (a) and (b) of Fig. 2 respectively. With p = and chiral limit, the two f's coincide and we have for the 
sum f(p, k) = 2c 2 ko(E — 77170). Comparing with Eq.(2.17) we get 

ImB h2 (E 2 ) = ^nidwiDldlwi (l, ~) (3.5) 

It can now be seen that both d and the Borel transforms of B 2 (E 2 ) and B^(E 2 ) are of 0(T A ). [The leading pieces 
from each of the two cuts are of 0(T 3 ), but they mutually cancel each other.] Thus the diagrams of Fig. 2 contribute 
neither to the pole nor to the cut to 0(T 2 ). 

The usual contribution of the irN intermediate state is given by the diagram of Fig. 3, for which f(p, k) — 0375^ — 
# + 777)75, c 3 = -(3/4)(A/^) 2 . For p = 0, the leading piece is given by / = 03(777 — JqE). Then the Borel 
transformed amplitude is 

A 2 T 2 e -™ 2 /A^ 2 

Finally the diagram of Fig. 4 gives a constant vertex correction to the nucleon pole, 

V6%E^^ + ^ ^ 

We are now in a position to write the complete spectral representation for the two point function to 0(T 2 ). 
Collecting results from this and earlier sections, the modified nucleon pole term is given by 
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E 2 
with 



(\T)2 

■ >( m + yoE), (3.8) 



32 F* 

while the Borel transform of the spectral representation due to irN branch cuts is, 

\2 rp2 -m 2 /M 2 

-{m{g\ 1) + {g\ + l) 7o £}_^^_ (3.10) 

The pole term immediately gives the results of Ref. 0: the effective nucleon mass does not shift to 0(T 2 ), while the 
shift in the residue is given by Eq.(3.9). The terms proportional to g\ in the branch cut contribution would have 
been absent, had the remainder R not been taken into account. The other remainder term S does not contribute to 
Q(T 2 ). 



IV. SUM RULES 



The other element needed in writing down the QCD sun rules is the short distance expansion of the product of 
nucleon currents. The quark content of the nucleon current, most suitable for the nucleon sum rules, is given for 
proton (i.e. i = 1, to be omitted henceforth) by |Tc|| , 

VD (x) = e abc (u aT (x)C ltl u»(x))( l5ltl d c (x)) D , (4.1) 

where C is the charge conjugation matrix. Here a, b, c are the color indices and D is a Dirac index. 

As already stated, we shall find the finite T sum rules, from which the vacuum parts are subtracted out, thereby 
eliminating the contribution of the unit operator. Since the leading thermal contribution of the gluon operators are 
of 0(T A ) and we are working to 0(T 2 ), the contributing operators are only uu, uu ilv t a uG^ va and uTiuuT2U up to 
dimension six. Of these the dimension five operator has zero coefficient in the OPE of the nucleon currents Q . Also 
the thermal expectation value of the four-quark operator turns out to be T- independent Thus it is only the 
operator uu which brings in any T-dependence on the operator side of the sum rules. We derive systematically all 
these results in Appendix B. 

Usually in writing sum rules in medium, one takes the parameters of the pole term as unknown, to be determined 
by the sum rules. Here these are A T and m T , the shifted pole position. Then the sum rules are obtained by equating 
the Borel transforms of the spectral representation and the operator product expansion. The 'subtracted' sum rules 
read, 

(A T ) 2 m T e-(™ T ) 2 /M 2 ^ + -^-(g 2 - \)— e -^IM^ ( IQ\ UU \Q\ K IL^ 

v ' M 2 16M 2 ' F 2 V 2 tt/ 8F 2 ' 

(A T)2 e -( mT )V^_L + J^tf A + l)^e- 2 /^= 0. (4.2) 

where the bar over the nucleon term denotes subtraction of its vacuum value. The factor K, given by 

.'2 



K = 1 



deviates from unity due to incorporation of the continuum contribution on the spectral side: the imaginary part of 
the leading operator contribution is assumed to saturate the spectral function from a threshold W onwards ^ . 
Let us evaluate the sum rules for m T and A T . Write 



i T 2 \ ( T 2 

m T = m{l + a—j , A T = A \l + b— 



where a and b satisfy, 



G 



M 2 2 32 vaA y 16m A 2 V 2tt 

-32^ + D. (4.4) 

The value of A may be obtained from the vacuum sum rules. Instead of using the numerical estimate |I, we may use 
one of the sum rules directly, which reads 

A 2 = _! f^l) (0\uu\0)Ke m2 / M2 (4.5) 
m \ 2tt J 

Then Eq.(4.3) becomes independent of the Borel mass M and one immediately gets, 

o = 0, b = -±(g 2 A + l), (4.6) 

reproducing the result of Ref. @ . 

It should be noticed that the prediction (4.6) is not as clean as it appears. Eq.(4.3) has M-dependence, though 
mild, within its allowed range of variation. The use of the vacuum sum rule for A 2 compensates this variation with 
its own to yield the above results. 



V. CONCLUDING REMARKS 



In this work we propose a method of constructing spectral representations for correlation functions in the medium. 
It is based on its Feynman diagrams which, to first order in the distribution function, consist of all 1-loop diagrams. 
This construction differs from the usual saturation scheme of taking only the single particle (with medium dependent 
mass and coupling) and two-particle intermediate states by certain 'remainder' terms contributed by the 1-particle 
reducible diagrams. 

We then use this spectral representation to write down the QCD sum rules for nucleon at finite temperature. They 
are shown to reproduce the temperature dependence of the nucleon mass and its coupling to the nucleon current, 
obtained earlier justifying simultaneously the spectral construction and the sum rules. 

These sum rules were obtained earlier by Koike Q , though his analysis of the spectral function is somewhat obscured 
by taking the chiral limit from the beginning. Of the two 'remainders', he calculates one explicitly and it so turns out 
that the other one, which he does not mention, is 0(T 4 ), when evaluated. As a result his sum rules remain unaltered 
by our method of construction of the spectral representation. But there are other sum rules, like the ones for the 
vector mesons p^Jl5| , which does have non-zero contributions from the 'remainder' terms, not taken into account so 

It will be observed that the Feynman diagram approach to constructing the spectral function automatically yields 
also the medium dependence of the mass and coupling of the single particle communicating with the current. It would 
thus appear that the only use of the sum rules in the medium is to rederive these results. However, the situation may 
not be so and one may well be able to extract new information on the matrix elements of operators. For example, 
consider the nucleon sum rules in nuclear medium. In this case there arises the nucleon matrix element of the four- 
quark operator, mentioned already in the introduction. Usually one relates it to the a-term by the approximation of 
ground state saturation. The sum rules, on the other hand, may give the value of this matrix element without such 
an approximation. 



APPENDIX A 



Here we derive a general formula for the discontinuity of 1-loop Feynman graphs encountered in this work. Each 
of these graphs consists of a nucleon and a pion propagator but with different vertices. Thus they are of the form 

Fu(p)=i J ^0 J /(p,fc)A 1 (fc)A 2 (p-fc) (A.l) 

where A12 are the 11-components of the corresponding propagator matrices [T^ ]; Ai(fc) is for pion, 
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A i(*0 = To o^— + ^n 1 (k )S(k 2 - f i 2 ), (A.2) 

and A 2 (p) is for nucleon, after extracting the spinor factor + m), 

A 2 (;p) = - ^— - - 27m 2 (;po)<5(p 2 - m 2 ) (A.3) 

p z — m z + te 

Here n\ and n 2 are the bosonic and fermionic distribution functions: ni(fco) = (e^ ' — n 2 (p ) = (e^ + 
The vertices and the spinor part of the nucleon propagator are all contained in f(p, k). 

Since we work at zero chemical potential and at T < /i, the nucleon distribution function n 2 is negligible compared 
to n\ for pion. However, we wish to retain both ri\ and n 2 to show the cancellation of their product in the final 
expression for the imaginary part. 

We find the imaginary part of Fn(p) by simply integrating out the time component of k^. For this purpose, we 
write the finite T propagators as 

A 1 (k) = {l + n 1 )D 1 (k)+n 1 D* 1 (k) (A.4) 
A 2 (p) = (1 - n 2 )D 2 {p) - n 2 D* 2 {p) (A.5) 



where 



k 2 — fi 2 +ie' p 2 — m 2 + ie 

The k integration over the integrands f(p, k)Di(k)D 2 (p — k) and f(p, k)Di{k)D 2 (p — k) are easily done by the residue 
theorem, getting simple poles in p on the real axis defined by the appropriate ±ie. Discarding the principal parts, 
we get the imaginary part as 

ImF n (p) 

/d 3 k 1 
— T - [{(1 + m)(l - n 2 ) - mn2}{/(wi)i(£! - wi - lu 2 ) + f(-u{)6{E + uj 2 )} 
(2tt) Aujiuj 2 

+{m(l - n 2 ) - (1 + ni )n 2 }{f{-^)5{E + Wl - c^ 2 ) + /K)<5(S - Wi + w 2 )}] (A.6) 

Here 



= \Jk 2 + v 2 , lu 2 = \J (p- k) 2 + m 2 . 



and rii = ni(wi) and n 2 = n 2 (w 2 ). For brevity, the argument of / shows only the value of the integrated variable fco- 
Observe the 'wrong' signs in the factors involving the distribution functions, because of which the product n\n 2 
does not cancel out. It is, however, possible to extract a factor of tanh(J3E/2) from each of the terms, by virtue of 
the associated delta functions, which 'corrects' the signs leading to their cancellation. Using a notation similar to 
Eq.(2.4) we get 

/d 3 k 1 
[(1 + m - n 2 ){/(wi)<5(£ - wi - oj 2 ) - f{-U!)S(E + wi + uj 2 )} 

+ K + n^ifi-u^SiE + lu 2 ) - f( Ul )5(E - lu, + lu 2 )}] (A.7) 

The positions of the branch cuts, where the imaginary parts are non- vanishing, are determined by the arguments of 
the (5-functions. 

We next consider the limit p = 0, in which we write the sum rules. The simplified cut structure in this limit are 
shown in Fig. 5. It is seen that the first and the third terms in Eq.(A.7) give rise to cuts for E > 0, and the second 
and the fourth terms for E < 0. The cuts in these two regions are related by symmetry under E — > —E. Restricting 
to cuts for E > and setting n 2 = we get 

/d 3 k 1 
-— 3 {(1 + m) f(u!)5{E - wi - wa) + m f(-e>i)8(E + wi - c^ 2 )} (A.8) 
(27r)° Auj\uj 2 

To evaluate the integral we have, for £7 > 0, 
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^S(E fwi- wa) = <H(£ T ^i) 2 - ^2} - ^(^i T w i) 



where wi = (E 2 

as given by the delta functions is always positive. We thus get 

ImF(E)\ E>0 = -L^gQa^Ql-^fiQi), 



ui + fi 2 )/2E. Note that u>\ is positive (negative) on the unitary (short) cut, but the value of ui 

(A.9) 



where g = 1 + n% iorE > m + fx and g — n\ for < E < m — ji. In writing the sum rules, we subtract the vacuum 
contribution (corresponding to the term '1' in <?). Then ImF(E) has the same expression on both the cuts. 



-(m+n) 


-(m- n) 


(m-n) 


/ E 

(m+Li) 







FIG. 5. The cut structure in the i?-plane 



APPENDIX B 



Here we derive systematically the known results of operator product expansion of nucleon currents using the 
co-ordinate space method [18-21 . We also collect the vacuum and the thermal matrix elements of the contributing 
operators 

The method consists in treating the gauge field as classical and Wick expanding the operator product. The basic 
element is then the contraction of two, say, u-quark fields, 



u(x) a A u(fi) a A> =~iS(x] 



aa 
AA< 



where the propagator S(x) satisfies 



ir{d^-iA^{x))S{x) =<5 4 (x), 

1 A" 



(B.l) 



in presence of the gauge fields A^(x) = gA"(x)^-, A™'s being the Gell-Mann matrices for SU(3) C . In the Fock- 
Schwinger gauge, defined by x^A^x) — 0, the gauge potential can be expanded in a series in the field strength 
and their covariant derivatives, 

Ap,{x) = — x a G afl (x) + • • • . 

The crucial step in this method is to solve Eq.(B.l) in a series at short distance |2(J 

S(x) = S (x) + S 2 (x) + ■ ■ ■ , 
where So 0*0 is the free propagator for massless quarks, 

1 ($AA> ,aa' 



S (x) a A a A> = 



2ir 2 (x 2 - ie) 2 
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and 52 (x) is the first non- leading piece at short distance, which is proportional to the gauge field, 

' l07T z (X z — It) 

These two terms suffice for our purpose. Note that Sq(x) is diagonal in color. 
With rj£,{x) given by Eq.(4.1), we have 

fj D ,(x) = (u a (xWC- 1 u(x) bT )(d c (x) 757v ) D >e abc . 

We are interested in the two point correlation function, 

T VD (x)fjD'(0) = ^^'"'^(C^AB^C-^^B^^DC^h'D' W^\(x,0) (B.2) 

where 

W$(x,0) = Tu(x) A u(x) b B d(x) c c u{0) a >{0) b B ,d{Q) c c, (B.3) 

is the operator product to be expanded into local operators with (singular) c-number coefficients as x^ — > 0. Its Wick 
expansion consists of three types of terms, 

w {a} _ j{a} jj{a} rrrM 
{A} _ 1 {A} + 11 {A} + 111 {A} 

corresponding respectively to one, two and three contractions of the quark fields. 

The single contraction terms may be obtained by contracting two rf's or two it's. In the latter case, there 
arises 4 terms, all of which are equal, as may be verified by interchanging the color indices and the Dirac indices in 
the terms and noting that C 7fl and 7 I /C~ 1 are symmetric matrices. Thus 

I { { a 2 } (x) = -iS(x)%' c , u{x)\u{x) B u{Q)i,u(Q%, - ±iS(x) a A a A , u{x) b B d{xf c u{Q) b Bl d{Q) c c'- (B.4) 
Similarly terms with two and three contractions give rise to 

Il\%{x) = ~2S(x) A a A ,S(x) bb BI d(0)c>d{x)° c - 4S(x) A a A>S(xy c c ' c , u{0)%u{x) b B (B.5) 

and 

IIl\ a A \{x) = -2iS(x) A a A ,S(x) bb B ,S(x)%'c> (B.6) 

Since we do not need pure gluon operators, we do not work with Eq.(B.6). The quark fields in Eq.(B.4) and (B.5) 
may now be expanded as 

u(x) = u(0) + xPDpuQS) + ^-D^D u u{fi) + ■ • ■ , (B.7) 

to get the series of local operators. Our task is now to project out the Lorentz scalar operators uu, 5 = ua^G^u 
and UT1UUT2U. We consider the vacuum matrix element for this purpose. 

Let us begin with lj^.(x). As we look for operators of dimension no higher than six, we replace S(x), u(x) and 
d{x) by the first terms in their expansions, 

l\% (x) = -iS Q (x)g c , u\ u b B u% u% - 4iS (x) A a A , u b B d c c u b B , & (B.8) 

where u\ = u A (0) etc. While projecting out the scalar part of the 4-quark operator, one simultaneously uses the 
approximation of vacuum saturation to relate it to the two-quark operator. For the first term we have 

(0\u A u b B ui,u b B ,\0) = (-S aa '6 bb 'S AA ,5 BB ,+S ab 'S ba 'S AB ,S BA ,) • (Q| ^J 0) , N = 12, (B.9) 

which gives rise to traces over 7-matrices and sums over e abc, s. It is simple to evaluate the traces after removing the 
C-matrix. We get 
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(0\Tr, D (x)fj D ,(0)\0) Itlst term = ^(0\uu\0) 2 ■ (B.10) 

A similar treatment to the second term results in a trace over an odd number of 7-matrices, giving 

(0\T r, D (x)fj D > (0)|0)j, 2nd term 

= 0. (B.ll) 

Next consider the twice contracted piece II j^,. The product of the two £(ai)'s may be expanded upto dimension 
two in the gauge fields to get 

C abc € a b C S(x)Xa'^( x )bB' — 2<5 CC S a aA'Sq BB' - S AA'^BB'^afl ~ ^0 BB'^A^A'Gap (B-12) 



In association with the matrices and7^C 1 in the first term of Il\ a \, the last two terms above are equal and we 



get its contribution to the operator product as, 



Tr) D (x)l] D ,(0) \n,lst term = -^{ibl^DC (TST'Oc 



D' 



tr(C^S ^C- l S^) ( dc>dc + x x d c >D x d c + ?-^d c ,D x D a d, 



c-Ux-Ucrac 



-tr{C ltl SoluC- l T af) T )d c ,G aP d c ] ■ (B.13) 
Note that the color indices in the operators are now summed. To project out the operator dd from Eq.(B.13), we use 

(0\d c >d c \0) = ^f-(0\dd\0), (B.14) 



to get 



2 

T T)D (x)f]D'{0) I O)//, 1st term, uu = 3j— g (0\dd\0)5DD< (B.15) 



ir^x 

We project out the other operator O5 by 



(0\d A D^D u d B \0) = — <j 5 B A1^ - ?Mba } (0\O 5 \0), (B.16) 



from which we also get 



(d A G au d B ) = ^-(a^) BA (0\O 5 \0). (B.17) 
It turns out that the two terms in Eq.(B.13) contributing to O5 mutually cancel each other, 

(Q\Tr) D (x)fj D ,(Q)\0) n , lst t erm,0 5 = 0. (B.18) 

We are yet to find the contribution of the second term of Eq.(B.5). On using Eq.(B.12), it can be seen to consist 
of terms containing the _B_B'-element of a matrix with an odd number of 7-matrices times operators like ub'Ub, 
UB'D\D a UB and u,B'G a 0UB- Also Eq.(B.14, B. 16-17) show that their projection on to the scalar operators produces 
a £>73'-element of an even number of 7-matrices, so we finally get a trace over an odd number of 7-matrices. We thus 
get § 

(0\Tr] D (x)fi D '(0)\0) IIi 2nd term, uu,Os =0 (B.19) 

Thus as far as the operators uu, O5 and the 4-quark operator are concerned, we have for the vacuum correlation 
function, 



i J d A xe l P x {0\T VD (x)fj D ,(0)\0) = P 2 H-p 2 ) (0|<ta|0)W - ~<0|<«|0) 



2 (B.20) 



Since the operators concerned are Lorentz scalars, one may think of getting the result for the thermal correlation func- 
tion by simply replacing the vacuum expectation values in the above equation by their ensemble averages. However, 
this would not be true for the four quark operator. For, in writing Eq.(B.9) we have already used the isospin (and 
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spin) structure of the vacuum matrix element, which differs from that for the ensemble average. A new calculation is 
thus necessary for the latter. 

To first order in the pion distribution function, the thermal average of an operator O is given by fll4| 



i=l,2,3 J v ; 



{k)\0\^{k)). 



(B.21) 



On using the soft pion methods, the pion matrix element may be reduced to the vacuum expectation value of a 
double commutator, 



Tr( P o) = (o\o\o)-±J ^1^1 E nUAQiom 



(B.22) 



where Q\ is the axial- vector charge, Q\ — J d 3 xA l (x). An elementary way to evaluate these commutators is to write 
Q\ in terms of quark fields, express the commutators in terms of anti-commutators of quark fields and replace the 
latter by their canonical values. They may then be vacuum saturated as before. 

We record first the well-known contribution of uu to the thermal nucleon correlation function, 



n(p) 



1 



T 2 



2tt J \~ 8F 2 
For the 4-quark operators in Eq.(B.8), we have 



tt) 1 - ^2 (OHO) p M~p 2 ) ■ i 



E (oi 

i=l,2,3 

(0\uu\0) 2 
N 2 



Qlu a A uW A'^ B 



|0) 



2{ZlAA>l { BB> + {l {, l5)AA,{l {> l 5 )BBiW a '5 bb ' - (A', a 



and 



E (0| \Ql\QluUc^B'd U C']} 10} 



i=l,2,3 



(0\uu\0) 2 
N 2 



-Hfr&B'-yCC ~ h°j5)BB'h a j 5 )cC'}5 bb '5 CC ' +i( 1 ° l5 )BC'(i a i 5 )cB'5 bc '5 cb ' 



(B.23) 



(B.24) 



(B.25) 



When these evaluations are inserted in the two point function, it turns out that the two operators bring contributions 
equal in magnitude but opposite in sign. As a result, there is no thermal contribution to 0(T 2 ) from the 4-quark 
operators, 



n ( ^)— "-rn , o • ^ j <o|««|o> 2 J* 



4 — quark 2 

3 



I 2 

F 2 



(B.26) 
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